Idealized glass transitions under pressure: dynamics versus thermodynamics 
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The interplay of slow dynamics and thermodynamic features of dense liquids is studied by exam- 
inining how the glass transition changes depending on the presence or absence of Lennard- Jones-like 
attractions. Quite different thermodynamic behavior leaves the dynamics unchanged, with impor- 
tant consequences for high-pressure experiments on glassy liquids. Numerical results are obtained 
within mode-coupling theory (MCT), but the qualitative features are argued to hold more generally. 
A simple square-well model can be used to explain generic features found in experiment. 

PACS numbers; 64.70.Q-, 62.50.-p 



The quest for identifying the physical mechanism be- 
hind the dynamical transformation of a liquid into an 
amorphous solid, the glass transition, has prompted 
many studies aiming to disentangle the dominant con- 
trol variables involved. It is recognized that the slow 
dynamics connected with the glass transition is univer- 
sal, but its connection to the underlying liquid structure 
is highly debated [H, [3, H, 0, IHl- Some argue in terms 
of a density effect called free or excluded volume; others 
attribute the main physics to energetic interactions and 
thermally activated processes. 

Experiments changing both temperature T and pres- 
sure P close to the transition are emerging to resolve such 
issues, but have brought contradictory results. Some 
find that temperature dominates glassy dynamics by far 
i, S i, H; some that it does not 0, [Hi, [H- 0th- 

equal influence 
some with temperature 
some with density g being more rele- 
3ll | . The difficulty of obtaining data over 
wide pressure ranges might be impeding: few studies, pi- 
oneered only in the 1990's [Ulli, go beyond 1 GPa. 
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Here we propose that to resolve the apparent contra- 
dictions, one needs to separate non-universal thermody- 
namic aspects, namely the equation of state (EOS) of the 
system, from universal dynamical features, viz. the slow 
relaxation. Specifically, we show how the presence or 
absence of attractive interactions affects the glass transi- 
tion, and how this emerges in different pairs of variables 
linked by the EOS: {g, T) (preferred by theory) vs. (P, T) 
(more amenable to experiments), yielding a transition 
that appears 'temperature-driven' in the latter. 

The Lennard-Jones (LJ) potential serves as a realistic 
interaction model: Vi,:i{r) = 4e[(r/(T)^^^ — (r/cr)~^], with 
dimensionless parameters g* = ga^ , T* = 1/ (/?e) , P* = 
Pa^ /e\ (3 = l/{kBT) with Boltzmann's constant k-Q. To 
study relative effects of entropy and energy, we compare 
the LJ glass transition with that of its purely repulsive 
(LJR) counterpart, li,jR(?') = Vi,j{r) -I- e > for r < 
2^/^(T, Vljr(?') = else. For the transition lines, mode- 
coupling theory (MCT) [H, HI] provides a reasonable 
qualitative description. Its transition temperature Tc is 



systematically above the experimental (calorimetric) one, 
Tg, but nevertheless serves as a good indicator for the 
change from high-T liquid-like dynamics to low-T glass- 
like one [3^ [36[. In MCT, Tc is the point where the 
glass form factor /(g), measuring an elastic response to 
the scattering spectrum, jumps discontinuously (usually 
from zero to a finite value) due to a bifurcation in [s^l 
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l-f{q) 2g4 J (27r) 

V{q,k) = S{k)S{p){q^ - qk/S{k) - qp/S{p)f/g^ with 
wave number q= |q|,p = |q — fc|. Equation ^ is solved 
numerically by an iteration scheme [37| . The interaction 
potential V{r) and temperature T enter only through the 
static structure factor S'(g), obtained in the hypernetted- 
chain (HNC) approximation [H, H^. S{q — 0) deter- 
mines the pressure through the EOS, 



f3P 



dgS{q = 0, g, T)' 



(2) 



integrated numerically for the LJ and LJR systems. We 
have checked that the quantitative error of the HNC clo- 
sure does not influence our results qualitatively. 

The MCT glass transitions for the two systems are 
almost identical, i.e., they depend very little on the pres- 
ence or absence of attractive interactions. As shown in 
Fig. [2 the lines in a T-g diagram coincide on the scale 
of the figure. In fact, most of the change in Tc{g) can 
be understood from the soft-sphere limit: for T ^ 00, 
^J(R)('') ~ n = 12, and the only control pa- 

rameter is an effective packing fraction, ipcS ~ (^/6)f?Ccffj 
where a^g = cr^ T* accounts for the soft core. The 
dash-dotted line in Fig. [T] corresponds to this soft-sphere 
glass transition, gc oc Tc where Eq. ^ yields ip^g « 
0.564. It clearly shows the same qualitative behavior 
as the Tc{g) lines for the full LJ system. Molecular- 
dynamics simulation data on the LJ glass transition col- 
lected in Ref. 0| and estimates from energy-landscape 
calculations of Ref. [4l| are reproduced in Fig. [T] as trian- 
gles. They agree with our results reasonably well, consid- 
ering the different definitions of the glass-transition point 
used in these works. 
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FIG. 1: Glass transition lines for the Lennard- Jones system 
with and without attractions: MCT resuhs for both coincide 
(sohd hne). Triangles are glass-transition points from sim- 
ulations jiol, |4l[. Dash-dotted line: soft-sphere asymptote, 
Tcie) ~ {^/&'Pch)'^q'^, Veil = 0.564. The dotted line indicates 
a gas-liquid spinodal. Inset: LJ (solid line), LJR (dashed), 
and square- well {S — 0.12; dash-dotted) potentials. 



Attractions have a crucial effect on the thermodynam- 
ics: they introduce a gas-liquid spinodal, whose estimate 
is shown in Fig. [T] (dotted line). The compressibility di- 
verges smoothly there, 1/S{q = 0) ^ 0. Hence, the 
resulting P{g,T) values in the hquid, Eq. are sig- 
nificantly lower with than without attractions, i.e., for 
T < Tcr and g > Qcr, where (pcr, ?cr) is the gas-liquid 
critical point. This even holds for approximations that 
fail to predict the spinodal (such as HNC), but replace it 
with branch points where S{q = 0) is large [i^. In other 
words, LJ-like attractions in high-density liquids simply 
affect the pressure: the contributions of all particles add 
up to a flat background that does not influence the forces 
[43], nor the (glassy) dynamics. 

Figure [2] demonstrates the marked effect the spinodal 
has on the glass transition line in the P-T diagram: 
although nearly indistinguishable in the T-g diagram. 
Fig. [l] the transition lines with (LJ) and without (LJR) 
attraction (thick solid and dashed lines) now appear qual- 
itatively different. The LJ log Pc-versus-log T line has a 
steep part not present in the LJR line at log^Q T* < 0.1, 
corresponding to T* < T*^. ^ lA. It stops at the spin- 
odal, restricting the glassy regime to densities g > Pcr 
(4^ . Simulation data from Refs. [i^, [3] (triangles in 
Fig. [2]) scatter between the LJ and LJR lines since the 
simulations used different truncations of the potential, 
yielding different equations of state. The data from 
Ref. |;40], where full {g, T, P) triplets are given for several 
truncations, collapse to a single curve in Fig. [T] within 
error bars, confirming our finding. 

The above results are easily understood within MCT: 
the glass transition in dense systems is driven by nearest- 
neighbor interactions (the cage effect), i.e., by features of 
the structure factor S{q) at q w lix ja. But the presence 
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FIG. 2: Lennard- Jones glass transitions with (thick line) and 
without (dashed) attractions in a pressure-vs-temperature 
representation, from MCT. Triangles: simulation data from 
Refs. pol . 13] (filled: LJ; open: purely repulsive; inverted: 
truncated LJ). Dash-dotted hne: LJ transition as P*ff = 
p,y-i/4^ Thin solid line: square-well model transition [5 = 
0.12, shifted as T* ^ l.ST*); dotted: square-well model spin- 
odal with critical point (cross) . Diamonds: experimental data 
for glycerol (Refs. [2^. l45j). see text. Arrows indicate isobars 
and isotherms used in Fig. [S] 



of attractions affects only S{q 0) [47|; a region strongly 
suppressed in the MCT integral, Eq. Merely the 

transformation p i— s- P is dominated by g — > effects. 

Since this transformation in general relies on straight- 
forward but cumbersome numerical calculations, we sim- 
plify matters by introducing a square- well (SW) potential 
as a 'cartoon' of the LJ model, showing first that the qual- 
itative features discussed above are still preserved. The 
SW model consists of a hard-sphere core and an attrac- 
tion of relative width 5: Vsw('') — — e for 1 < r/a < 1+5. 
Here, a mean-spherical approximation (MSA) for S{q) 
'is'] allows to integrate Eq. ([2]) analytically, greatly sim- 
plifying calculations. The MCT hne for 5 0.12 [33] is 
shown in Fig. [21 rescaled as T* l.ST* to account for 
the different S{q) approximation that mainly induces a 
shift in the T scale 
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We identify two generic regimes for both the LJ and the 
square- well lines: T* ^ 00 is the hard-sphere limit where 
the glass transition occurs along an isochore Pc ^ k-gT] 
this limit is approached for T* > 1 (regime I). It is also 
present in the LJ system, provided one corrects for soft- 
core effects, P*ff = p*j^-3/n^ ^j^g dash-dotted hne 
shows. For 0.1 < T* < 1 the steep log P*-versus-logT* 
line discussed above is found in both models (regime II). 
Its position along the T* axis scales with the gas-liquid 
critical temperature (T^*. ~ 5 for the MSA-SW). For 
T* ^ 0, the SW model shows a low-density regime cor- 
responding to densities g < pcr (a cross marks the point 
where gc — Per), where the dynamics itself strongly de- 
pends on the potential depth e. This attraction-driven 
regime may be connected with colloidal gelation and 
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FIG. 3: Typical relaxation times r as functions of density: 
MOT results for the SW system {S = 0.12; bottom and 
left axes), along the isobar P* — 0.316 (solid line) and the 
isotherm T* = 0.3 (dashed), marked by thick arrows in Fig. [2] 
Symbols (upper /right axes) are viscosity data for isopropy- 
Ibenzene [13] for the isobar P = 0.1 MPa (filled) and the 
isotherm T = 293K (open). Inset: MCT resuhs for P* = 8.91 
(solid) and T* = 0.775 (dashed); thin arrows in Fig.[2l 



is absent in the Lennard- Jones model and in common 
molecular glasses. Relevant for typical glass formers at 
MPa pressures is regime II, as pointed out in a recent 
study [i^: experiments reveal steep P{Tg) curves that 
are incompatible with the hard-sphere-like regime I. Di- 
amond symbols in Fig. [2] exemplify this for glycerol: ex- 
perimental Tg{P) data from Refs. [2S 
according to e/fce = 500 K and e/cr^ 



45[ was mapped 



2.5 GPa, just to 
demonstrate qualitative agreement (and absorbing the 
quantitative difference between Tg and Tc in the map- 
ping). 

Our discussion of glass-transition lines has direct im- 
plications for dynamical quantities, as the latter are ex- 
pected to depend strongly on the distance to the transi- 
tion. For example, a 'thermodynamic scaling' has been 
observed for many glass formers, that involves gT~^ as a 
scaling variable [8|, ISll . [53,[5l|,[53 and can be interpreted 
as an effective density-dependent interaction energy 0. 
7 is an empirical, effective exponent: even in the LJ sys- 
tem, 7 ^ 1/4 = 3/n [5^. In agreement with this, we 
find that the LJ transition line can be well fitted for all 
T* < 3 by g*{T) oc T'^, where 7 = 0.15 . . .0.23. Partly, 
this merely mirrors that effective power laws can be used 
to fit the potential in the respective V{r)/€ ~ T* range. 
If attractions are present, the effective 7 also depends on 
their details, as strikingly demonstrated by the SW sys- 
tem, where we do not find the 7 — » expected from the 
n —^ 00 hard-sphere limit. 

One can characterize the relative effects of tempera- 
ture and pressure on the glass transition by monitoring 
the viscosity 77 along isotherms and isobars, using density 
as a parameter. Experiments [1, [13, Ei [13, El El HE, 
H, [ii H m, Ei, [23, [H, ^ usually find stronger varia- 



tion along an isobar (varying T) than along an isotherm 
(varying P). This is a natural consequence of our sce- 
nario, evidenced in Fig.[3]by the viscosity rj* in SW units, 
calculated from the standard Green-Kubo expression in 
MCT approximation [s^. Consider regime II: changing 
T along the P* = 0.316 isobar (solid line) corresponds 
to a more direct approach to the glass transition line as 
compared to changing P along the T* = 0.3 isotherm 
(dashed) ; cf. the thick arrows in Fig. O Clearly, -q* di- 
verges over a smaller density interval along the isobar. 
The agreement with experiment is semi-quantitative, as 
the symbols, reproduced from Ref. [17[, demonstrate. 
Only in regime I (inset of Fig. [Sj thin arrows in Fig. [2|) 
do P and T excert equal influence. Not all transition 
points are equal: MCT predicts nonuniversal amplitudes 
and shapes for relaxation spectra that change along the 
transition line, invalidating a 'temperature-pressure' su- 
perposition principle. But the changes are small enough 
to make it appear to work, explaining why some exper- 

some with restric- 

MM. 
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iments find it 

Mm 



some not at all 



tions ^ ^ ^ ^ 

Temperature- and density-effects are often quantified 
by a ratio of activation energy and enthalpy [o, [HI, , 
Ev/Hp ^ ((91ogT/ar-i)v/(91ogr/ar-i)p, trivially 
related to E = [d g / dT)r{dT / d g) p . Large \E\ signify a 
temperature-driven transition, small \E\ a density-driven 
one. E consists of a glass-transition part, and a purely 
EOS-driven one. According to Fig.[Tl the isokinetic term 
r = {dg/dT)r does not depend on LJ-like attraction; the 
EOS term t = {dT/dg)p however changes. Estimating 
the latter through the SW-MSA expression, we find that 
in fact, t decreases in the vicinity of the spinodal: the 
measure \E\ increases with P, indicating a growing in- 
fiuence of temperature at higher pressures. Such trends 
have been found in experiments and argued to be at odds 
with the expectation that the transition becomes hard- 
sphere like at high P [i^. According to our model, they 
are dominantly thermodynamic effects. 

A similar conclusion holds for the pressure depen- 
dence of 'fragility' often used to classify how quickly re- 
laxation times diverge. Recent work debated its rela- 
tion to g — > quantities such as elastic constants, the 
above energy-enthalw ratio or the effective exponent 7 
0, [3, [3, [3, H [H [53. If true, MCT predicts a pres- 
sure dependent fragility that is a nonuniversal feature of 
S{q 0), less so of the glassy dynamics [g^. This may 
also hint towards why the presence of hydrogen bonds 
drastically changes the high-pressure behaviour regard- 
ing fragility and thermodynamic scaling (sij l . 

In conclusion, MCT glass-transition lines for the 
Lennard- Jones system and for the same system truncated 
to be purely repulsive are nearly identical in a g-T plot. 
Fig- HI Yet, they appear quite different in a P-T diagram. 
Fig. [2] The difference can be understood as unrelated to 
the glass-transition mechanism itself, but to a difference 
in thermodynamic behavior only. If one accepts that the 
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glass transition is a primarily dynamic phenomenon, it 
will not be altered by sufficiently long-ranged, LJ-like, 
attractions. However, the equation of state, determining 
the pressure of the system, will change. The Pj,-versus-T 
curve hence has two regimes if attractions are present. In 
the very-high pressure regime I, it is essentially a density- 
driven fluid-glass transition: Pc oc T, after correcting 
for soft-core effects. In regime II, identified as the ex- 
perimentally relevant one, the existence of a gas-liquid 
spinodal leads to Tc{P) curves with a much weaker P- 
dependence: this could be called a "temperature-driven" 
liquid-glass transition; but "temperature driven" is not 
equivalent to "attraction dominated" . Key qualitative 
features can be understood with the help of the square- 
well system as a better tractable model. 

Discussing the glass transition in terms of "tempera- 
ture vs. pressure" might obscure the physics responsible 
for it, focusing too much on different thermodynamics of 
different glass formers. For example, a change in compo- 
sition in metallic glass formers greatly changes the ther- 
modynamics, but has little effect on the slow dynamics 
. Experiments probing the pressure range P ^ 1 GPa 
seem desirable to test the picture proposed here. 
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